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NERVA DYNAMIC ANALYSIS METHODOLOGY 
-SPRVIB- 


I. INTRODUCTION 

This report documents the general dynamic computer code called SPRVIB 
(Spring Vib) developed at Westinghouse Astronuclear Laboratory in support of 
the NERVA (Nuclear Engine for Rocket Vehicle Application) Program. Using normal 
mode techniques, the program computes kinematical responses of a structure 
caused by various combinations of harmonic and elliptic forcing functions or 
base excitations. Provision is made for a graphical type of force or base 
excitation input to the structure. 

II. SUMMARY AND CONCLUSIONS 

The usefulness of a dynamic computer code can be determined by examining 
four main features : 

1) the method of solution 

2) the type and generality of the forcing functions used 

3) the flexibility of the input/output routines 

4) size limitations of the program caused by storage and/or computational 
limits 

The computer code, SPRVIB, should then be examined for all of these features. 

A description of the required input format and a listing of the program are 
presented, along with several examples illustrating the use of the program, in 
Reference 1. SPRVIB is written in Fortran IV computer language for use on the 
CDC 6600 or the IBM 360/75 computers. 

SPRVIB was used by ANSC for the dynamic analysis of the NERVA 400D and 
400D' engine configurations. For the 400E configuration it was not used at all 
due to the loss of in-house computing facilities, the limited size of the program, 
and the adoption of the NASTRAN program as the standard analysis tool for the 
structural dynamics group. 


D. F. Vronay 
26 May 1972 


N8120R: 72-030 
Page 2 


III. TECHNICAL DISCUSSION 

SPRVIB uses a classical normal mode solution in which the individual re- 
sponse is computed as a sum of all the modal responses. This type of solution 
has two main advantages: 

1) It completely describes the structure and often enables one to esti- 
mate the frequency dependence of the structure. This, of course, is not possible 
with a numerical solution. 

2) It enables the exact solution to a linear problem to be determined 
in closed form if the forcing functions can be described by elementary functions. 
In these cases, roundoff type errors do not accumulate. 

A normal mode solution also has two disadvantages: 

1) At each time step in a transient solution to be printed out, conver- 
sion must be made from the normal coordinates to the physical coordinates of the 
problem. 

2) Only proportional damping or modal damping type solutions are tract- 
able for general programming. 

The first disadvantage does not appear to be limiting in these days of the 
high speed computer. If the forcing functions are easily describable, such as 
harmonic or elliptic functions, the solution obtained is exact at each time step 
and the step size can be determined by convenience or plotting limitations. On 
the other hand, the time step size of the ordinary numerical type solutions de- 
termines the accuracy and stability limitations of the solution. Furthermore, 
with the possible exception of physical damping instruments, the damping con- 
stants used in dynamic analysis are usually crude estimates of the actual damping 
As such, it would appear at least as correct to use estimates of modal damping 
constants which often may be obtained from test results. Hence, the apparent 
limitations are often not applicable to practical problems. 

SPRVIB offers the following wide variety of possible forcing functions: 

1) Free vibrations with arbitrary initial displacement and velocities 
at each generalized mass. 
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2) Arbitrary harmonic forces F Q sinwt acting at each or all of the gener- 
alized mass points. Forces with different amplitudes and frequencies can be 
simultaneously applied at the various points. 

3) Elliptic type forcing F^expCBt) can act at each or all of the points 
similar to that described above. The parameters F q and 8 are real parameters 
which may have different values in as many as four consecutive time regions. 

Using this type of representation, a shock or explosive type of forcing can be 
easily described with sufficient accuracy to obtain quite accurate estimates of 
even velocities and accelerations. In addition, the harmonic and elliptic loads 
may be superimposed at any point. 

4) Base harmonic excitation of the supports. 

5) Point force input at any mass point. In this case the force is repre 
sented as a set of data points. 

6,7,8) Base motion in which the support acceleration (velocity, displace- 
ment) is described by a set of data points. 

The effectiveness of the input and output features of any program greatly 
affect the time a user must spend in setting up a problem and the practicality of 
its use. A limiting feature of most dynamic computer codes available at this 
time is the requirement that the element matrices, such as the stiffness matrix, 
be generated by the direct stiffness method or some other similar method. SPRVIB 
will generate the element matrices for discrete mass elements connected by rigid 
bars, linear and rotational springs, and viscous damping elements. Besides the 
regular input options, provision is made to input only the non-zero symmetric 
elements of sparse matrices. This provision was used exclusively for the NERVA 
dynamic analysis. This option is very useful for large matrices. All the input 
can be in a punched card format or some matrices can be read in from magnetic 
tape. 

An option is available to store all the matrices and results from the 
topological portion of the program, including the mode shapes and frequencies, 
on tape so that this section of the program will not have to be redone at a later 
date for further load cases on the same problem. Any transient information which 
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can be printed out can also be stored on tape for possible future plotting. Be- 
sides the regular steady state response, provision is made to output the results 
of each harmonic or elliptic load acting separately as well as the combined effect 
for each load case. In this way the transient effect of each load may be deter- 
mined with a coarse time step while the effects of the combined loadings can be 
determined with a finer time step. One, two, or all of each of the displacements, 
velocities and accelerations may be output at each time step. 

At the present time SPRVIB is limited to 90 degrees of freedom. The 
storage requirement of the program with the 90 degrees of freedom is 175,000 o . 

O 

The detailed program write-up, including a listing of the program, is 
contained in Reference 1. The remainder of this report illustrates a typical , 
application of the program to the 400D' engine analysis. 

Figure 1 shows the dynamic model of the 400D' engine, consisting of mass- 
less beams connecting a series of concentrated masses so distributed as to approxi 
mate the mass distribution of the engine. Table 1 is a list of the degrees of 
freedom and what physical interface they represent. N 

The input to SPRVIB consisted of the diagonal mass and damping matrix, the 
stiffness matrix and a specification of the input/output control parameters. 

The output consisted of an echo of the mass, damping and stiffness matrices 
a list of the non-zero entries by vector locations, and a row by row printout of 
the final stiffness matrix. This is followed by the calculated output consisting 
of the damping constants, and finally the eigenvectors and corresponding natural 
frequencies. These modes are then used, together with the transformed input, to 
perform, in this example, the transient response analysis due to a base input 
acceleration of 0.3g. Two types of output are presented. First, a summary of the 
peak response at each degree of freedom. This is then followed by the actual time 
history of the coordinate displacements, velocities, and accelerations. 

The output from this sample run appears as Appendix A. Only a portion of 
the time history response is shown. 

IV. REFERENCES 

1. "SPRVIB A General Purpose Dynamic Computer Code Using Normal Mode 
Solution", Westinghouse Astronuclear Laboratory Report WANL-TME-1940 , June 1971 
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TABLE 1 

• - 




DEGREES OF FREEDOM 

Point 

X_ 

Y_ 

0_ 


A 

1 

2 

3 

Stage 

B 

4 

5 

6 

Interface (Stage/Upper Thrust Structure) 

C 

7 

8 

9 

Upper Thrust Structure/Gimbal 

D 

10 

11 

12 

Gimbal Pivot Forward 

E 

13 

- 

- 

Actuator //I Forward End 

F 

14 

- 

- 

Actuator //I Aft End 

G 

15 

- 

- 

Actuator #2 Forward End 

H 

16 

- 

- 

Actuator #2 Aft End 

I 

17 

18 

19 

Gimbal Pivot Aft 

J 

20 

21 

22 

Gimbal/Lower Thrust Structure 

K 

. 23 

24 

25 

Lower Thrust Structure - 

L 

26 

27 

28 

Lower Thrust Structure 

M 

29 

30 

31 

TPA #1 

N 

32 

33 

34 

TP A #2 

0 

35 

36 

37 

Shield 

P 

38 

39 

40 

Lower Thrust Structure 

Q 

41 

42 

43 

Lower Thrust Structure/Pressure Vessel 

R 

44 

45 

46 

Pressure Vessel 

S 

47 

48 

49 

Pressure Vessel 

T 

50 

51 

52 

Pressure Vessel 

U 

53 

54 

55 

Pressure Vessel 

V 

56 

57 

58 - 

. Pressure Vessel/Nozzle/Nuclear Subsystem 

CD 

59 

60 

61 

Nuclear Subsystem 

AB 

62 

63 

64 

Nuclear Subsystem 

W 

65 

66 

67 

Nuclear Subsystem 

AY 

- 

68 

- 

Stage 

EE 

69 

70 

71 

Nozzle Throat 

FF 

72 

73 

74 

Nozzle Divergent 

GG 

75 

76 

77 

Nozzle Torus/Extension 

HH 

78 

79 

80 

Nozzle Extension 

II 

81 

82 

83 

Nozzle Extension 

JJ 

84 

85 

86 

Nozzle Extension 
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DEGREES OF FREEDOM (Cont.) 


Point 

X_ 

Y 

0_ 


KK 

87 

88 

89 

Nozzle Extension/Destruct Subsystem 

*LL 

90 

91 

92 

Destruct Subsystem 



*The Destruct Subsystem parameters are set to zero to eliminate it. 
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TYPICAL RUN 


SAMPLE OUTPUT 
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